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Optomechanical response of a non-linear mechanical resonator 
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We investigate theoretically in detail the non-linear effects in the response of an optical/microwave 
cavity coupled to a Duffing mechanical resonator. The cavity is driven by a laser at a red or blue 
mechanical subband, and a probe laser measures the reflection close to the cavity resonance. Under 
these conditions, we find that the cavity exhibits optomechanically induced reflection (OMIR) or 
absorption (OMIA) and investigate the optomechanical response in the limit of non-linear driving of 
the mechanics. Similar to linear mechanical drive, an overcoupled cavity the red-sideband drive may 
lead to both OMIA and OMIR depending on the strength of the drive, whereas the blue-sideband 
drive only leads to OMIR. The dynamics of the phase of the mechanical resonator leads to the 
difference between the shapes of the response of the cavity and the amplitude response of the driven 
Duffing oscillator, for example, at weak red-sideband drive the OMIA dip has no inflection point. 

We also verify that mechanical non-linearities beyond Duffing model have little effect on the size of 
the OMIA dip though they affect the width of the dip. 


I. INTRODUCTION 

In optomechanics, light and mechanical motion are 
coupled together by the radiation pressure of photons 
trapped in an optical cavity. Optomechanics has the po¬ 
tential to provide access to the quantum limit of mechan¬ 
ical motion, and in recent years, has been a field that has 
been undergoing rapid development^. It already resulted 
in a number of ground-breaking experiments, including, 
for example, observation of radiation pressure shot noise^ 
or optomechanical squeezing of light^’^. Following the 
first observation of quantum nature of nanomechanical 
resonator®, experiments in microwave optomechanical® 
and cavity optomecanica/’® architectures demonstrated 
that mechanical resonators can be brought to the quan¬ 
tum regime. The next step will be to achieve reliable 
quantum manipulation and to use mechanical resonators 
as quantum memory elements and quantum transducers. 

An important tool in optomechanics is optomechani¬ 
cally induced transparency (OMIT), an effect analogous 
to electromagnetically induced transparency in quantum 
optics®’^®. In an OMIT experiment, the cavity is strongly 
driven by a drive (pump) laser at the red-sideband — 
the driving frequency is red-shifted from the cavity res¬ 
onance by the frequency of the mechanical resonator — 
and the transmission through the cavity is measured by 
a probe laser close to the cavity resonance. Due to the 
interference, the transmission exhibits a narrow peak ex¬ 
actly at the cavity resonance. The width of the peak is 
determined by the mechanical relaxation rate, and the 
OMIT peak is typically much more narrow than the cav¬ 
ity resonance. The observation of the OMIT peak serves 
as the signature of optomechanical coupling and can be 
used to qualitatively determine the coupling in the exper¬ 
iment. OMIT was first observed in optomehanical exper¬ 
iments in Ref. 12. It has been heavily used for charac¬ 
terization of optomechanical systems, in particular, for 
extraction of the value of optomechanical couplingme- 
chanical resonators and optomechanical coupling, in both 
opticaF’^®d3 microwave^^’^® realizations. In single¬ 


port cavities, like the one used in Ref. 15, one can only 
measure reflection. A peak/dip in the reflection coeffi¬ 
cient is referred to as optomechanically induced reflec¬ 
tion/absorption (OMIR/OMIA) 

In most optomechanical experiments, both cavity and 
mechanical resonator were linear. However, non-linear ef¬ 
fects are conceptually important. In the classical regime, 
a driven non-linear oscillator exhibits bistable behavior, 
which may strongly affect the properties of the system. 
Moreover, in the quantum regime, non-linear effects are 
essential for creation of non-classical states of a non-linear 
resonator. There are three sources of non-linear behavior 
in cavity and microwave optomechanics. First, the ra¬ 
diation pressure interaction is inherently non-linear, but 
typically a cavity is strongly driven by a laser, and in this 
case the interaction can be linearized^. There were theo¬ 
retical proposals to use the radiation pressure interaction 
for quantum manipulation of mechanical resonators^®, 
which require strong coupling regime of operation. Sec¬ 
ond, a cavity can be made non-linear. Whereas it is diffi¬ 
cult to achieve for optical cavities, microwave cavities can 
be made non-linear by adding Josephson junctions. Self- 
sustained oscillations caused by non-linear Lorentz force 
backaction^^ and cavity-induced Duffing behavior of a 
mechanical resonator^® were demonstrated in Josephson- 
based devices the dc regime. Third, mechanical res¬ 
onators are inherently non-linear. This is best manifest 
in low-dimensional resonators such as carbon nanotubes 
and graphene flakes^®. In this Article, we concentrate on 
non-linear effects in mechanical system. 

Non-linear mechanical effects in optomechanical archi¬ 
tecture were experimentally studied in Ref. 20 (see also 
Ref. 21). These articles investigated OMIT in an un¬ 
dercoupled optical cavity and discovered that the OMIT 
feature is still in place, and the shape of the OMIT peak 
repeats the shape of the frequency response of the driven 
non-linear (Duffing) oscillator and exhibits bistability re¬ 
sulting in a hysteretic behavior. 

Recently, we performed experimental studies of non¬ 
linear effects in OMIA®® in the same overcoupled mi¬ 
crowave cavity as Ref. 15. The results we found have a 
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FIG. 1. Schematic of an optomechanical cavity consisting of 
a static (left) and a movable (right) mirrors. The movable 
mirror acts as a mechanical resonator. It is coupled to the 
cavity through the radiation pressure. 


qualitatively different lineshape from that of Ref. 20 — 
the shape of the OMIA dip in the reflection coefficient of 
the cavity was distorted by non-linear effects, however, 
the shape is different from the response of a driven Duff¬ 
ing oscillator, in particular, it does not have an inflection 
point. 

In this Article, we provide a theoretical analysis of non¬ 
linear effects on OMIR/OMIA feature and explain the 
difference between Refs. 20 and 22. We demonstrate that 
whereas the OMIR/OMIA feature persists in all situa¬ 
tions and exhibits the bistability at strong enough probe 
powers, its shape is different depending on how the cavity 
is coupled to an external circuit. We show that the non- 
Duffing shape of the OMIA response is due to the effect 
of the phase of the mechanical resonator imprinted on 
the microwave field. We also analyze the OMIR/OMIA 
for a blue-detuned drive. 

The paper is organized as follows. In Section II, 
we develop a general theory of non-linear effects in 
OMIR/OMIA based on input-output relations — a stan¬ 
dard technique in quantum optics, and derive the OMIA 
for a red-sideband driven overcoupled cavity. In Section 
III we generalize the results to the case when the drive 
is detuned from the red sideband, and in Sec. IV we 
demonstrate what happens for an undercoupled cavity, 
connecting our results with Ref. 20, and for a blue- 
sideband drive. In Section V, we consider non-linear ef¬ 
fects beyond the Duffing appoximation; The quadratic 
term in the force and the non-linear dissipation. Section 
VI presents the conclusions. 


II. MODEL OF A DRIVEN NON-LINEAR 
CAVITY 


We consider a single-port cavity, which is coupled via 
the radiation pressure to a non-linear resonator as shown 
in Fig. 1. The Hamiltonian of the system is 


H = huica^d + 




2mef / 2 4 

+hGxd^d -I- ihy/rfK [sm(t)a^ — s*^(t)d] . 


( 1 ) 


Here a^(a) is the creation (annihilation) operator for the 
cavity mode with the cavity frequency Wg. The resonator 
with the effective mass mg// and the resonance frequency 


flm is described by the operators of mechanical displace¬ 
ment X and momentum p operators. The mechanical 
non-linearity is described by the Duffing term with the 
strength a. The radiation pressure term contains the 
optomechanical coupling constant G = dtOc/dx between 
the mechanical and cavity modes, which corresponds to 
the shift of the cavity frequency due to the motion of 
the resonator. The last term describes the coupling of 
the input signal to the cavity. The cavity coupling pa¬ 
rameter is given by 77 = k^/^kq + Kg) = /Cg/re, where 
kqj Ke and K denote the intrinsic, the external, and the 
total dissipation rates, respectively. The input signal 
Sin{t) = -I- consists of the strong drive 

field with the normalized field amplitude Sd and drive fre¬ 
quency bJd and weak probe field with the normalized am¬ 
plitude Sp and probe frequency LOp. The amplitudes Sd 
and Sp are square roots of the corresponding drive/probe 
power divided by fvjjc and have dimensions of S^ 

and 5^ have the meaning of the incident photon flux 
measured in photons per second. In the remainder of the 
article we will refer to the amplitudes Sd and Sp simply 
as drive and probe helds. 

Writing the Heisenberg equations of motion in the 
rotating-wave frame of the drive frequency^’^^ and adding 
dissipation of the cavity and mechanical modes, we de¬ 
rive the quantum Langevin equations for our system as 
follows, 

^a(t) = d{t) - iGx{t)d{t) + y/rjKSin{t)e'‘‘^’^*,{2) 

= —tiGd\t)d{t) — ax^{t), (3) 


where A = Wg — is the detuning between the cavity 
field and drive field and T^ is the damping rate of the 
resonator. 

In the steady-state, disregarding the probe field, the 
time derivatives vanish, and the static solutions for the 
intra-cavity field and mechanical displacement obey the 
following algebraic equations, 


Vw a 

a = —= - Sd, 

-iA + K/2 

meffO^x + ax^ + hGa^ = 0 , 


( 4 ) 

( 5 ) 


where A = A — Gx is an effective cavity detuning includ¬ 
ing the frequency shift due the static mechanical displace¬ 
ment. 

Since the probe field is much weaker than the drive 
field, following the standard methods of quantum optics, 
we can rewrite each Heisenberg operator as the sum of its 
steady-state mean value and a small fluctuation, which 
has zero mean value 


d{t) = a + 5a{t) and 5x{t) = x + 5x{t). ( 6 ) 

In this case, the steady-state values are governed by 
the drive power and the small fluctuations by the probe 
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power. Then, keeping only the linear terms of the fluctu¬ 
ation operators in the radiation pressure term (disregard¬ 
ing 6a^{t)6a{t) and 6a{t)Sx{t)), we obtain the linearized 
quantum Langevin equations 

^5a{t) = (lA — ^)Sa{t) — iGadx{t) + ^JrjKSpe~''^*', (7) 

meff + ^m^Sx{t) + Q'^5x{t)^ = -hGa{Sa{t) 

+Sa\t)) — a{Sx^{t) + Sx'^6x{t) + SxSx‘^(t)), (8) 

with Q, = ujp — ujd being the detuning of the probe field 
from the drive. In order to solve this system of equations 
we introduce the following Ansatz: Sait) = -I- 

We are interested here in the resolved sideband 
regime {k <C flm) and close to the blue (A = or 
red (A = —ilm) sideband of the drive field, meaning 
the lower sideband A'^ is far off-resonance and can be 
neglected. Upon substituting this Ansatz into Eq. (7) 
we derive the amplitude A of the cavity field 

A= -=-^-A_!— £ -. (9) 

-f(A + U) + f -i(A + fi) + § 


Next, we substitute this into Eq. (8). This is a non¬ 
linear equation which can be solved in the same way as 
one solves the forced Duffing oscillator^^. We rewrite the 
equation in terms of the dimensionless time r = Qmt and 
frequency ui = 


-^Sx{t) +u}'^Sx{t) = -(w^ - 1)(5 x(t) - t^-^Sx{t) 


dr^ 


rirn dr 


hGd 

a 


{Sx^{t) + 3x^Sx{t) + SxSx'^{t)). 


( 10 ) 


The terms proportional to x^ and xSx'^{t) are dropped 
since they create a very small static force. To fa¬ 
cilitate analysis of this non-linear equation, the small 
asymptotic parameter £ <g; 1 is introduced. Then, as¬ 
suming weak non-linearity, weak damping, weak forc¬ 
ing with frequency close to the natural frequency of 
the mechanical resonator, one can perform transforma¬ 
tions: (x/uiefy SCx/uief^m/^m ^ 

hGdA/mf.f —>• ehGdA/mef a-nd (w^ — 1) —>■ 
£(w^ — 1). This indicates that all parameters on the r.h.s. 
of eq. (10) are small, and we can construct the solution 
as a power series using the method of multiple scales 


5x{t) = xo{to,ti) + exi{to,ti) + C>(e^), (H) 


with the fast timescale to = t and the slow timescale 
ti = £T. Applying two timing leads to the transformation 
of the first and second time derivative with respect to 
time r 


jL 

dr 

dr^ 


A A 

dto ^ dti 


— Dq + eDi, 


D^ + 2eDiDo+0{e^). 


( 12 ) 

(13) 


Using relations for the derivatives in the equation of mo¬ 
tion for mechanical resonator and collecting orders 0(1) 
and 0(£) yield the pair of differential equations 


D^xo +oj‘^xo = 0, 


(14) 




DqXi + UJ Xi = w — 1 -I- -= 


h 




Xo 


—2DiDoXo — — — DqXq — 


—*(A -|- ^ 

hGdy/riKSpe~^‘^*°/ruef 

— i(A -|- UlVlm) -\- ^) 

■(a:o + 3x^X0). (15) 


^e// 




The first equation is a simple harmonic oscillator with the 
general solution xo{toAi) = -|-A*(ti)e®‘^‘‘’. 

Substituting xo to the r.h.s. of eq. (15) and collecting 
terms proportional to the we obtain the force that 

drives the l.h.s. at its resonance frequency. In order for 
the perturbation correction xi to not diverge, the sum 
of all collected terms must be zero. It gives rise to the 
condition for the slowly varying amplitude X(ti), 


2iujX ■ 


-1 


— z(A + UjQ^rn) A ^ 


hGay/^Sp/meff^l, 

—A. — -^- 

—2(A -t- LU^rfi) ^ 


-3- 


a 




■{XX* +x'^)X = 0. 


A 


(16) 


This is our main result, which in the remainder of the 
Article we will simplify looking specifically at the blue or 
the red sideband of the cavity field. 


A. Red sideband 


For the red sideband, the drive field is fixed close the 
lower motional sideband A = — 0^ + <5, where <5 is a 
small shift (detuning) from the mechanical frequency. 
The probe field is slightly detuned from the cavity reso¬ 
nance frequency, D = -I- A' — <5, where A' is a small 

detuning. We also define = hG'^d^/2meffilm and 
F = hGdy/r[KSp/meffilm- Then Eq. (16) is approxi¬ 
mated as follows 


2iX + 2 



\ 


X 


F 


-iA' 


- 3- 


meff 


Hr 


-{XX* + x^}X = 0. 


(17) 


Introducing a polar form of the complex amplitude 
Al(ti) = and separating real and imagi¬ 

nary part {2iX = ib+b(f>) gives a system of two first-order 
differential equations for the real amplitude b and phase 
(j). These equations describe a slow time evolution of 
the amplitude and phase. To find the equilibrium points 
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of this slow flow we look at the solution with constant 
amplitude(& = 0) and phase((^ = 0) 


FA'sin((^) — —Fcos{4)) = g'^bA' 


+ A'' 


3a 


‘IvTief 


-(6V4 + x2 )-A' + c^') &, (18) 
/ 


F — sin(((i) + FA' cos{(j)) 


h. 


(19) 


The results for the amplitude of the mechanical res¬ 
onator is shown on Fig. 2 at zero detuning of the drive 
{5 = 0) and a > 0 as a function of the probe frequency 
for two probe powers. We see that the behavior is exactly 
the same as for the Duffing oscillator^^. At low driving 
(low probe power) the response of the amplitude is the 
same as for a linear oscillator and results in a Lorentzian 
peak. If the probe power increases, the resonant curve 
becomes asymmetric and develops instability, shown as 
hysteresis. 



FIG. 2. (Color online) The amplitude of the mechanical res¬ 
onator rescaled by the probe field Sp = 10^s“^^^(red) and 
Sp = 3 X 10®s“^/^(blue) as a function of the frequency A'/k 
driven exactly at the red sideband. For the non-linear fre¬ 
quency response (blue) the hysteresis occurs. The blue arrow 
up shows the direction of the jump in the amplitude for de¬ 
creasing frequency and the blue arrow down shows the direc¬ 
tion of the jump for increasing frequency. Most parameters 
are taken from the experimental paper by Singh et. al}^ 
with the cavity coupling parameter g = 0.777, except for 
the intra-cavity number of photons = 7.6 x 10®, mechan¬ 
ical damping r,„/27r = 200 Hz, and non-linearity stregth 
a = 2 X 10^"^ kg . 


The amplitude of the cavity field can now be found 
using real amplitude and phase of the mechanical oscil¬ 
lator. 


-iGabe *'^/2 

-f(A + D) + f ^ -i(A + 0) + f ■ 


( 20 ) 


B. Blue sideband 


For the blue sideband, we follow the same steps as 
for the red sideband by taking the frequencies D = 
uip — ujd = —-b A' — i5 and A = ujd — uic — Gx = 
+ b. Consequently, the equations for the ampli¬ 
tude (18) and phase (19) can be adapted by chang¬ 
ing sign in the term on the r.h.s. of Eq. (18) by 

- -f A'2) bA' ^ -h -b A'2) bA' and in Eq. (19) 
by (4 -b A'2) r^5/2 ^ - (4 + A'") T^b/2. 


C. Reflection coefficient 


In order to study OMIA or OMIR, we need to look at 
the cavity reflection of the probe field. It is defined by 
the ratio of the output and input field amplitudes at the 
probe frequency. The output field of a single port cavity 
can be found using the input-output relationship 


'5oui(t) — -b . (^1) 


Hence, the reflection coefficient of the probe field is given 

by 


l^iil 



( 22 ) 


For a linear mechanical resonator, the amplitude A is 
linear with the probe field Sp, because the fluctuations 
of the mechanical resonator 6x are linear with Sp. Hence, 
the reflection coefficient is independent of the probe field. 
For a non-linear resonator the relation is more complex. 

At zero detuning (5 = 0 and for a linear resonator, the 
shape of OMIT as the function of the detuning between 
drive and probe exactly repeats the shape of the response 
of the amplitude of the mechanical resonator. We see 
from Eqs. (18), (19) that this is not the case for a non¬ 
linear resonator. The main reason is that OMIT is not 
only determined by the dynamics of the amplitude b, but 
also by the dynamics of the phase shift (j), which in turn 
depends on the frequency. This dependence is essentially 
an effect of back-action of the cavity on the resonator. 
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FIG. 3. (Color online) Optomechanically induced absorp¬ 
tion window found from the cavity reflection coefficient ISnl, 
which corresponds to the mechanical amplitude shown in Fig. 

2 . 

The non-trivial shape of reflection is iiiustrated in Fig. 
3 for the same parameters as Fig. 2. Note that in the 
remainder of the Articie we wiii continue to use the same 
parameters to visuaiize the anaiysis uniess mentioned 
otherwise. Oniy stabie parts of the curve are shown; the 
bistabiiity resuits in jumps between different branches 
of the reflection. It is seen that the reflection close to 
the cavity resonance is suppressed, indicating OMIA for 
a single-port cavity. The shape of the OMIA dip for 
stronger probe (Duffing oscillator) is bistable, however, 
it does not correspond to the shape of the resonator re¬ 
sponse (Fig. 2), for instance, it does not have an inflec¬ 
tion point. 

Fig. 4 illustrates the dependence of the OMIA dip on 
the probe field Sp for the same parameters as in Fig. 2. 
The red curve corresponds to the field Sp = 
when the dynamics of the mechanical resonator is lin¬ 
ear, and is the same as the red curve in Fig. 3. The 
blue curves, from left to right, correspond to increasing 
probe power, and the resonator shows non-linear behav¬ 
ior. We see that upon increasing the probe frequency, the 
width of the OMIA dip increases and the depth slightly 
decreases with the probe power. It does not correspond 
to the Duffing oscillator feature, where the peak height 
stays unaffected by the drive power. 

In the rest of the Article, we investigate the effect of 
different parameters at the response of the cavity. 


III. OMIA DETUNING 

In this Section, we look at the effect of the detuning of 
the drive S. The results are shown in Fig. 5 for the same 
parameters as Fig. 2 and various values of S. 

Fig. 5 a and b give the general picture of the cav¬ 
ity response in a broad frequency window for positive 
(a) and negative (b) detunings at weak probe field, 
Sp = 10^5“^/^. The detuning breaks the symmetry of 
the response with respect to the center of the cavity res¬ 
onance even in the linear case, shifting the OMIA dip. 


FIG. 4. (Color online) OMIA response for different probe 
fields. The first red curve shows linear response at the probe 
field Sp = The dashed grey line represents the peak 

height of the linear response, in order to compare it with other 
probe fields. The blue curves start on the left from the probe 
field Sp = 1 X and equally increase the probe field 

by Sp = 0.5 X step up to the last curve on the right, 

which corresponds to S'p = 6 x 


We see indeed that the drive detuning shifts the narrow 
OMIA dip away from the minimum of the broad cavity 
resonance. This can be also seen from Eq. (18): Indeed, 
S enters there in the combination S + 3ab'^ / {Snif, ff^m)- 
We have seen already that increasing the probe power 
(increasing h) shifts the OMIA minimum to the right, 
therefore positive 5 must shift it to the right as well, and 
the negative detuning must shift the peak to the left. 

Fig. 5 c and d show the structure of the dip for the 
same parameters as a and b, respectively. The shape is 
deary asymmetric and is different from the Lorentzian 
shown in Fig. 3. The flank which is closer to the cen¬ 
ter of the cavity resonance (the right/left one for neg¬ 
ative/positive detuning d) is sharper than the opposite 
flank. The OMIA dip at the detuning i5 is a mirror re¬ 
flection of the OMIA dip at the opposite detuning —5. 

Fig. 5 e, f and g, h show the cavity response at 
a stronger probe fields, Sp = 3 x and Sp = 

5 X when the non-linear nature of the response 

is well pronounced. First, OMIA dips at 6 and —6 
are no more mirror images of each other. For the for¬ 
ward sweep in e and f, if the non-linear term is positive, 
a > 0, the OMIA dip at negative detuning —6 is less 
pronounced than the dip for positive detuning S. On 
the contrary, for the backward sweep in g and h, the in¬ 
crease in the probe power results in the decrease of the 
depth for both, positive and negative, detunings. For a 
linear system, we would expect a symmetry of the sig¬ 
nal 5'ii(a;p — ujc'jS) = Sii{—ujp + Wc;—<5), which would 
make left and right panels of Fig. 5 mirror reflections of 
each other. In the presence of non-linearity, the curves 
are pulled in the same direction. Thus, the symmetry is 
broken, and therefore the OMIA dip is pronounced for 
the negative detuning. 
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Linear response 




(c) 




Non-linear response 






FIG. 5. OMIA at non-zero detuning. All panels on the 
left show OMIA for the negative detuning <5 = —0.2k; all 
panels on the right are for the same parameters but the pos¬ 
itive detuning 5 = 0.2k. (a), (b) The full cavity reflection 

coefficient in the broad frequency range for the low probe 
field Sp — when the cavity can be considered to 

be linear, (c), (d) Zoom of the linear OMIA peak, (e), (f) 
Non-linear OMIA window of the forward frequency sweep for 
the two probe fields: Sp = 3 x (dashed blue line) 

and Sp = 5 X (solid green line), (g), (h) Non-linear 

OMIA of the backward frequency sweep for the same probe 
fields as in (e) and (f). 


IV. NON-LINEAR RESPONSE MAP 


So far, we concentrated on red-sideband drive and over¬ 
coupled cavity, p > 1 / 2 , and a moderate drive power. 
In this situation there is an OMIA dip which is again 
shown in the left curve of the left bottom cell of Fig. 
6 . If the drive power increases, the cavity response Sn 
instead of a minimum shows a maximum. This is an 
OMIR peak, as demonstrated in the right panel of the 
same cell of Fig. ( 6 . Similarly for OMIA, at weak probe 
powers the OMIR peak is Lorentzian, and if we increase 
the probe power, the peak shifts to the right and broad¬ 
ens. Note that whereas the shape of the OMIA dip for 
weaker drive is different from the response of the me¬ 
chanical resonator (no inflection point), the shape of the 
OMIR peak for stronger drive repeats the shape of the 
resonator response. 

In the bottom right cell of Fig. 6 we show the results 
for the same overcoupled cavity driven at the blue side¬ 
band. We see that the OMIR peak develops at any drive 
power. For low probe powers the peak is Lorentzian, and 
for higher probe powers it shifts to the left and broad¬ 
ens. The shape of the peak is related to the shape of the 
response of the mechanical resonator, Fig. 2. 

The qualitative difference between red- and blue- 
shifted drive for OMIR/OMIA has already been ana¬ 
lyzed in Ref. 15 for the linear dynamics of the mechan¬ 
ical resonator. Indeed, the minimum value of the re¬ 
flection coefficient is a non-monotonous function of the 
coupling parameter p: It decreases with rj for an un¬ 
dercoupled cavity, reaches zero for an optimally coupled 
cavity 77 = 1 / 2 , and increases with rj for an overcou- 
pled cavity. The red/blue-sideband drive modifies the 
cavity linewidth kq such that it becomes kq ± dp^/Fm, 
where the upper/lower sign corresponds to the red/blue- 
detuned drive. Thus, for an overcoupled cavity increasing 
the drive intensity (proportional to the number of pho¬ 
tons in the cavity) for the red-sideband drive takes the 
cavity towards the undercoupled limit, and the behavior 
changes dramatically when the effective coupling crosses 
the point 77 = 1/2, crossing over from OMIA to OMIR. 
In contrast, the blue-sideband drive only takes the cav¬ 
ity to even stronger overcoupled regime, and there are 
no qualitative changes. For an undercoupled cavity, the 
roles of red- and blue-sideband driving are swapped. 

Now we turn to the non-linear dynamics. From Eq. 
(18) we see that the signs of 6 ^ and A' are opposite for the 
red-sideband drive and the same for the blue-sideband 
drive. This means that since the OMIR peak shifts to the 
right from the cavity resonance for the red sideband, it 
shifts to the left for the blue sideband, in full accordance 
with Fig. 6 . 

The top row of Fig. 6 shows OMIA/OMIT for an un¬ 
dercoupled cavity, 77 < 1/2. The red-sideband drive (top 
left corner) was previously studied in the experiments^®. 
For any drive powers, there is an OMIR peak, and the 
shape of the peak correponds to the shape of the Duff¬ 
ing resonator response. The top right corner shows the 
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FIG. 6 . Non-linear response map schematically shows the position and the shape of the OMIA/OMIR features for the 
undercoupled (77 = 0.223), in the first row, and overcoupled (77 = 0.777), in the second row, cavity. The two columns rep¬ 
resent red and blue sideband. When the intra-cavity number of photons rid = ( in the plots (a), (b), (d), and (f) ) the 

reflection coefficient is shown for the probe field Sp = with the solid(red) line, 5p = 3 x with the dashed(blue) 

line and Sp = 5 x with the dash-dotted(green) line. For the larger intra-cavity number of photons rid = 25a^ (in the 

plots (c) and (e) ) there is a switch between OMIA and OMIR , then the reflection coefficient is plotted for the following probe 
fields Sp = with the solid(red) line, Sp = 3 x with the dashed(blue) line and Sp = 5 x with the 

dash-dotted(green) line. Each plot corresponds to the different value of the Duffing non-linearity stregth a = 2 (a,d), 0.5 (b,f), 
0.08 (c,e) xlO^'^ kg 


cavity response for the blue-sideband drive. It is sim¬ 
ilar to what happens in an overcoupled cavity for the 
red-sideband drive: At low probe powers, one has an 
OMIA dip with the line shape different from the Duffing 
response of the resonator (no inflection point), whereas 
for stronger driving, an OMIR peak develops with the 
Dufhng-like shape. The position of the peak is shifted to 
the left of the cavity resonance with increasing the probe 
power. 

V. BEYOND DUFFING 

To complete our results, we consider two more non¬ 
linear effects not included in the Duffing model (1) — the 
quadratic term in the displacement and the non-linear 
damping. Both can be included in the treatment of Sec¬ 
tion II. 

One phenomenon is the additional term in the 


energy of the mechanical resonator. It produces the force 
proportional to the coordinate x breaking the symmetry 
of the oscillator potential and shifting the equilibrium 
point. Here, we restrict ourselves to the case /3 > 0. 
Fig. 7 shows the results for the overcouples cavity driven 
exactly at the red-sideband, ^ = 0. The curve for /3 = 0 
corresponds to the OMIA dip shown in Fig. 3. We see 
that the only effect of finite /3 is to shift the position of 
the OMIA dip to the left of the cavity resonance — in 
the direction opposite to the shift due to increasing probe 
power. The shape and the depth of the dip are almost 
unaffected. 


Next, we consider the effect of non-linear damping, 
adding the term —fix^x to the equations of motion for 
the resonator. Without the non-linear term, a = 0, this 
would be the van der Pol resonator^^, obeying the equa- 



























































FIG. 7. Evaluating the contribution of to the en¬ 
ergy of the mechanical resonator on the OMIA dip. The 
quadratic non-linearity = 0 is shown by solid(red) line, 
/? = 3 X 10*^^ kg by dashed(blue) line, and /3 = 

6 X 10^^ kg by dash-dotted (green) line for the fixed 

probe field S'p = 3 x The l3x^ term induces a force- 

dependent shift of the equilibrium frequency, but otherwise 
does not change the non-linear dynamics or optomechanical 
response. 



FIG. 8. (Color online) OMIA response including nega¬ 
tive non-linear damping for different probe field. The red 
curve shows the linear response with zero negative non-linear 
damping for the probe field Sp = The 

dashed grey line represents the peak height of the linear re¬ 
sponse. The next blue curves include negative non-linear 
damping g — 3x They start with the probe field 

Sp = 1 X and end with Sp = 4 x 10®having 

an equally spaced step in the probe field of 0.5 x 10®s“^^®. 
The effect of the negative non-linear damping term consid¬ 
ered here is to increase the depth of the OMIA dip in the 
optomechanical response with increasing driving force on the 
mechanical resonator. 


tion of motion 



where F is an external force. For sufficiently large oscil¬ 
lation amplitudes the system would become unstable far 
from the equilibrium point > F^), however, here we 


are interested in the situation when the non-linear term 
in the damping is a small correction. 

Fig. 8 shows the effect of the non-linear damping on 
the OMIA dip in the same situation as Fig. 3. There is 
an enhancement of the dip for the increasing probe power 
and positive g. Correspondingly, for negative g the dip 
becomes deeper with increasing probe power. 

VI. CONCLUSIONS 

In this Article, we systematically investigated the ef¬ 
fect of non-linearities in a mechanical resonator on the 
OMIR/OMIA feature of a one-port optomechanical cav¬ 
ity. The main result is summarized in Fig. 6. We see 
that depending on whether the cavity is overcoupled or 
undercoupled, whether it is red- or blue-sideband driven, 
and whether the drive power is sufficiently strong, the 
non-linear effects either result in OMIR with the shape 
repeating that of a driven Duffing oscillator (an analog 
of Ref. 20), or OMIA with a different shape, as observed 
in Ref. 22. The higher is the probe power, the wider 
is the feature and the more pronounced is the effect of 
non-linearity, including the width of the hysteresis range. 
The depth of the OMIA feature is not significally affected. 
Furthermore, we investigated the effect of other factors 
on the OMIA feature and found that the detuning of 
the drive shifts the OMIA dip and makes it significantly 
asymmetric, the quadratic term in the oscillator force 
shifts the position of the dip without considerably affect¬ 
ing its shape, whereas the non-linear dissipation enhances 
the dip. These conclusions will help to interpret the re¬ 
sults of the cavity response in optomechanical cavities 
with a non-linear mechanical component. 

Let us now discuss the limitations of our theoretical 
model. First, our theory is valid in the weak coupling 
regime, g <C k. Whereas this seems to be sufficient for ex¬ 
isting experiments with non-linear resonators^^’^^, multi¬ 
photon strong coupling has already been achieved in the 
linear mechanical limit We linearize radiation-pressure 
interaction, and thus we can also describe multi-photon 
strong coupling regime using a similar framework. Gen¬ 
eralization to the single-photon strong coupling regime is 
however difficult. Furthermore, we only considered the 
situation when the probe power is much weaker than 
the drive power. It would be relatively easy to con¬ 
sider a complementary situation when the probe is much 
stronger than the drive, however, we have use the method 
of scale separation, and it fails when the drive and the 
probe power are comparable. Equations still could be 
solved numerically in this case. 
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